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ABSTRACT 

We investigate the possibility that near future observations of ultra-high-energy cosmic rays (UHECRs) can 
unveil their local source distribution, which reflects the observed local structures if their origins are astrophysical 
objects. In order to discuss this possibility, we calculate the arrival distribution of UHE protons taking into account 
their propagation process in intergalactic space i.e. energy losses and deflections by extragalactic magnetic field 
(EGMF). For a realistic simulation, we construct and adopt a model of a structured EGMF and UHECR source 
distribution, which reproduce the local structures actually observed around the Milky Way. The arrival distribution 
is compared statistically to their source distribution using correlation coefficient. We specially find that UHECRs 
above 10 19,8 eV are best indicators to decipher their source distribution within 100 Mpc, and detection of about 
500 events on all the sky allows us to unveil the local structure of UHE universe for plausible EGMF strength 
and the source number density. This number of events can be detected by five years observation by Pierre Auger 
Observatory. 

Subject headings: cosmic rays — methods: numerical — IGM: magnetic fields — galaxies: general — 
large-scale structure of the universe 

1. INTRODUCTION 

The origin of ultra-high-energy cosmic rays (UHECRs) 
above 10 I9 eV is one of the most intriguing problems in as- 
troparticle physics. Akeno Giant Air Shower Array (AGASA) 
found statistically significant small-scale cluste rings of ob- 
serve d UHECR events with large-scale isotropy dTakeda et al.l 
fl999h . The AGASA data set of 57 events above 4 x 10 19 eV 
contains four doublets and one triplet within separation an- 
gle of 2°. 5, consistent with the experimental angular resolu- 
tion. The chance probability of observing su ch multiplets un- 
der a n isotropic distribution is only about 1 % ([Havashida et al.l 
2000). A combination of the results of many UHECR ex- 
periments (including AGASA) also revealed eight doublets 
and two triplets within 4° on a totally 92 events above 4 x 
10 19 eV (lUchihori et al.1 120001) . These multiplets suggest that 
the origins of UHECRs are point-like sources. For identi- 
fication of UHECR sources, arrival directions of UHECRs 
have been observed in detail by High Resolution Fly's Eye 
(HiRes) and Pierre Auger Observatory (Auger). However, so 
far, these experiments have reported no sign ificant clustering on 
the arrival distribution above 4 x 10 19 eV dAbbasi et al.l l2005a: 
iMollerach etai]|2007l) . 

Recently, several classes of astrophysical objects in many 
literature has tested for positional correlations with observed 
arrival directions of UHECRs. The correlations with BL Lac 
objects were discussed on the assumptions of smaller deflec- 
tion angles of UHECRs than the experimental angular re solu- 
tion and/or neutral primaries (Tinvakov & Tkachev 2001), and 
in consid eration with the deflection b y Galactic magnetic field 
(GMF ) dTinvakov & Tkachevl 12002b . iGorbunov & Troitskvl 
(120051) considered various classes of powerful extragalactic 
sources for the correlation study includin g small corrections 
of UHECR arrival directions by GMF. lHague et al.1 d2007l) 
discussed the correlation with nearby active galactic nuclei 
(AGNs) from RXTE catalog of AGNs. However, these stud- 
ies have not taken into account UHECR propagation in ex- 
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tragalactic space. UHECRs above 8 x 10 19 eV lose a signifi- 
cant fraction of their energies by photopion production in col- 
lision with the cosmic micro wave background (CMB) pho- 
tons during their propagat ion dBerezinskv & Grigorieval ll988; 
lYoshida & Teshimal Il993l) . Thus, UHECRs have horizons, 
which are the maximum distances of their sources that UHE- 
CRs can reach the Earth, even if their energies are below 
8 x 10 19 eV at the Earth. The positional correlations between 
arrival directions of UHECRs and their so urce candida t es out - 
side the horizons are not significant. (In Hag ue etail ((2007), 
only nearby AGNs within the horizons are considered.) In ad- 
dition to the UHECR horizons, deflections due to extragalactic 
magnetic field (EGMF) are also important since extragalactic 
cosmic rays are propagated for a much greater distance than 
in Galactic space. Propagation process of UHECRs should be 
c onsidered in such correlatio n studies. 

lYoshiguchi et al.l (|2003a) investigated the correlation be- 
tween the arrival distribution of UHECRs and their source dis- 
tribution taking into account UHECR propagation in intergalac- 
tic space with a uniform turbulent EGMF whose strength is 1 
nG and coherent length is 1 Mpc. The authors adopted a source 
distribution with 10~ 6 Mpc -3 that reproduced the local struc- 
tures and the AGASA results. They concluded that detection 
of a few thousand events above 4 x 10 I9 eV reveal observable 
correlation with the sources within 100 Mpc. 

However, a uniform turbulent field is not realistic EGMF 
model. Faraday rotation measurements indicate magnetic 
field strengths at the /iG level within inner reg ion (~ cen- 
tral Mpc) of galaxy clusters (iRronbergl |1994|) . The evi- 
dence for synchrotron emis sion in numerous galaxy clusters 
( Giovannini & Ferettil |2000[) and in a few cases of filaments 
( Ki m et al.ll 19891: iBagchi et al.ll2.002b also seems to suggest the 
presence of magnetic fields with 0.1-1.0/iG at cosmologi- 
cal structures. Several numerical simulations of large-scale 
structure formation have confirme d these magnetic structures 
dSigl et alJfcOOHlDolag et al1l2005l) . 
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Based on these studies, in recent years, we calculated prop- 
agation of UHE protons in a structured EGMF which well 
reproduces the local structures actually observed and sim- 
ulated their arrival distributions with several normalizations 
of EGM F strength and severa l number density of UHECR 
sources (Takami &Sat3|2007|). We constrained the source 
number density to best reproduce the AGASA results. As 
a result, 10~ 5 Mpc -3 is the most appropriate number den- 
sity, which is weakly dependent on EGMF strength. (In rec- 
tilin ear propagation, s imilar number density is also obtaine d 
in iBlasi & de Marcol d2005l) : iKachelriess & Semikoj d2005l) ) 
However, this has large uncertainty due to the small number 
of observed events at present. 10~ 4 Mpc~ 3 and 10~ 6 Mpc~ 3 are 
also statistically allowed. Therefore, it is useful to deliberate 
the correlation between the arrival distribution and the source 
distribution in the case of these number densities. Note that 
we revealed in the paper that this uncertainty will be solved by 
future increase of detected events. 

In this study, we calculate the arrival distribution of UHE- 
CRs, taking their propagation process into account, and investi- 
gate the correlation the arrival distribution and their source dis- 
tribution in the future. A structured EGMF model and source 
distribution which can reproduce the local universe actually ob- 
served are adopted. The source number density and the EGMF 
strength are treated as parameters since these have some uncer- 
tainty. A goal of this study is that we understand the number of 
observed events to start to observe the UHECR source distribu- 
tion by UHECRs and how much the correlation is expected in 
the future. 

Auger has already detecte d more events abo ve 10 19 eV than 
those observed by AGASA dRoth et al .1120071) . Nevertheless, 
the event clustering has not observed, as mentioned above. It 
might be due to EGMF and/or GMF strong enough not to gen- 
erate the multiplets or statistical fluctuation for small number 
of observed events at highest energies. In any case, we should 
predict and discuss how the arrival distribution reflects UHECR 
source distribution. 

Chemical composition of UHECRs is very important for the 
correlation. If UHECRs are heavier components, magnetic de- 
flections are larger and the correlation is worse. One of observ- 
ables for study of UHECR composition is the depth of shower 
maximum, X max , which can be measured by fluorescence de- 
tectors. Its average value (X max ) is dependent on UHECR com- 
position and energy. HiRes reported that composition of cos- 
mic rays above 10 19 e V is dominated by p rotons as a result 
of X mH measurement ( Abb asi et al.l l2005bi) . Recent result by 
Auger is com patible to the H iRes result within systematic un- 
certainties dUnger et al.ll2007l) . However, they concluded that 
the interpretation of (X max ) distribution is ambiguous because 
of the uncertainties of hadronic interaction at highest energies. 
Thus, UHECR composition at highest energies is controversial 
at present. Despite that, in this study, we assume that all UHE- 
CRs are protons since composition above 10 19 eV has proton- 
like feature. 

This paper is organized as follows: in section [2] we pro- 
vide our models of UHECR source distribution and a struc- 
tured EGMF. In section [3] we explain our calculation method 
for the arrival distribution with UHECR propagation and sta- 
tistical method. In section |4] The results of the correlation be- 
tween the arrival distribution of UHE protons and their source 
distribution. We summarizes this study in section|5] 



2. SOURCE DISTRIBUTION AND MAGNETIC FIELD 

In this section, our models of UHECR source distribution 
and a structured EGMF are briefly explained. These models are 
almost the same as those i n our previous work . More detailed 
explanations are written in Taka mi et al. (2001). 

These models are constructed from Infrared Astronomical 
Satellite Point Sour ce Catalogue Redshift Survey (IRAS PSCz) 
catalog of galaxies dSaunders et al.l[2000i) . This is a catalog of 
flux-limited galaxy survey with large sky coverage (~ 84% of 
all the sky). Thus, this is appropriate galaxy catalog for the pur- 
pose. The selection effects for observatio n are corrected wit h 
luminosity function of the IRAS catalog dTakeuchi et al.l l2~003). 
We use a set of galaxies after the correction within 100 Mpc 
(called our sample galaxies below) for construction of our mod- 
els since small number of galaxies can be observed above 100 
Mpc. We adopt Vt m = 0.3, tt A = 0.7 and H Q = 71 km s" 1 Mpc" 1 
as the cosmological parameters. 

We assume that subsets of our sample galaxies with specific 
number densities are UHECR source distributions. 10~ 4 , 10~ 5 
and 10~ 6 Mpc -3 is considered as the source number density. 
For the source number densities, we randomly select galaxies 
from our sample galaxies with probabilities proportional to the 
absolute luminosity of each galaxy. This method allows con- 
structing source distributions to reflect the cosmic structures. 
Source distribution above 100 Mpc is assumed to be isotropic 
and luminosity distribution of galaxies follows the luminosity 
function. We take the sources until 1 Gpc into account. All 
sources are assumed to be have the same power for injection of 
UHE protons. 

Our EGMF model also reproduces the local structures ac- 
tually observed around the Milky way. Several simulations 
of cosmological structure formation with magnetic field have 
found that EGMF roughly traces ba ryon density distribution. 
dSigl etal.ll200l iDolaeet al.l l2005h . Our structured EGMF 
model is constructed with simple assumptions in addition to 
these results. We constructed matter density distribution from 
our sample galaxies with a spatial resolution of IMpc equal 
to the correlation length of our EGMF model, l c . In each 
cell, strength of the EGMF is related to matter density, p, as 
|-Begmf| oc p 2 / 3 and a turbulent magnetic field with the Kol- 
mogorov spectrum is assumed. The strength of EGMF is nor- 
malized to B = 0.0, 0.1 and 0.4/iG at a cell that contains the 
center of the Virgo Cluster since EGMF strength has some un- 
certainty as mentioned in section 1 . These three normalizations 
are used for investigation of the spatial correlation between the 
arrival distribution and the source distribution. Since the EGMF 
model and the source distribution are constructed from the same 
galaxy sample, UHECR sources are in the magnetized struc- 
ture. Note that volume of about 95% has no magnetic field in 
our EGMF model. 

3. METHOD OF CALCULATION 

3.1. Calculation of arrival distribution of UHE protons 

Once a source distribution is given, the arrival distribution 
of UHE protons can be calculated by calculating propagation 
of protons from their sources to the Earth. In this section, we 
describe fundamental processes on proton propagation in inter- 
galactic space and our calculation method of the arrival distri- 
bution. 

UHE protons lose their energies in collision with 
CMB during their propagation in intergalactic space 
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dBerezinskv & Grigorieval 119881: lYoshida & Teshimal Il993l) . 
Higher energy photon backgrounds (e.g. infrared, optical, ul- 
traviolet etc.) are neglected in this study since we treat protons 
above 10 19,6 eV. Such protons remain nearly unaffected by the 
higher energy background photons because of their relatively 
small number. Protons above 8 x 10 19 eV lose their energies by 
photopion production, p + 7 — ► tt+X. This reaction has large 
inelasticity (~ 30%) and relatively short energy-loss length 
of about a few ten Mpc at z = 0. Protons with such energies 
cannot reach the Earth from distant sources. Thus, the pho- 
topion production predicts sharply spectral steepening around 
8 x 10 19 eV, so-called Greisen-Zatsepin-Kuz'min ( GZK) steep- 
ening <[Greisi3[l966t ^m^^&^^^^\l96^. The GZK 
effect is essential in considering the correlation with distribu- 
tion of relatively distant sources. We adopt the energy-loss 
length which is calculated by simulatin g the photopion pro duc- 
tion with the event generator SOPHIA (Mu cks et alj|2000h . On 
the other hand, protons below 8 x 10 19 eV lose their energies 
mainly due to Bethe-Heitler pair creation, p + j — ► p + e + + e~. 
This has small inel asticity (~ 10~ 3 ). We adop t the analytical 



fit function given bv lChodorowski et alj d 19921) to calculate the 
energy-loss rate on isotropic photons. An adiabatic energy loss 
due to the expanding universe also makes protons lose their 
energies. The energy-loss rate is written as 



-^=-~E=-H [Q m (l +z) 3 + n A ] 



1/2, 



(1) 



These three energy-loss processes are treated as continuous pro- 
cesses in our calculation. 

Trajectories of protons are deflected by EGMF. Protons in- 
jected from their sources to the Earth cannot be led into the 
Earth straightforward. It wastes much CPU time to calculate the 
propagation of protons which cannot arrive at the Earth in order 
to construct the arrival distribution. For solving such problem, 
we suggested a new calculation method of the arrival distribu- 
tion which is an ap plication of the backtracking propagation 
dTakami et alj|2006b . In this method, protons with the charge 
of -1 are ejected from the Earth and then we calculate their tra- 
jectories in intergalactic space taking into account magnetic de- 
flections and energy-loss(gain) processes. Such trajectories are 
regarded as those of protons from extragalactic space. We cal- 
culate only trajectories of protons which can reach the Earth. 

In order to simulate the arrival distribution, 2,500,000 pro- 
tons (with charges of -1) with dN /d\og 10 E cx Const, from 10 19 
to 10 21 eV. are ejected isotropically from the Earth. We cal- 
culate their trajectories until their propagation time exceeds 
the lifetime of Universe or their energies reach 10 22 eV, which 
corresponds to the maximum acceleration energy at UHECR 
sources. For each source distribution, we calculate a factor for 
a trajectory of j th particle, which represents the relative prob- 
ability that the j th proton reaches the Earth, 

P fFfwV 1 dN/dE g (d i , j ,Eg i ) dE g 

Here i labels sources on each trajectory, while Zij and dij, are 
redshift and distance of such sources. E„' is the energy of a 
proton at the i th source, which has the energy E at the Earth. 
Thus, dN /dEg(d i j,E g i ) cx E g ~ 2 ' 6 is the energy spectrum of UHE 
protons ejected from a source whose distance is dij. dE g /dE 
represents correction factor for the variation of the shape of the 
energy spectrum through the propagation. 

We randomly select trajectories according to these relative 
probabilities, P se i e c, so that the number of the selected trajecto- 
ries is equal to the required event number. The mapping of the 



ejected direction of each particle at the Earth is the arrival distri- 
bution of UHE protons. If we have to select the same trajectory 
more than once to adjust the number of selected trajectories, 
we generate new events whose arrival direction is calculated 
by adding a normally distributed deviation with zero mean and 
variance equal to the experimental resolution to the original ar- 
rival direction. The experimental resolution is assumed to be 
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3.2. Statistical Quantities 



In order to investigate statistically the similarity between the 
arrival distribution of UHECRs f e and the source distribution 
f s , a correlation coefficient between the two distributions is de- 
fined as 

~ {feJs) ^(^l— ( 3 ) 



where 



P (f a j b )=j2 



VP<JeJe)P<Js,fs) 

fa(j,k)-f a f b U,k)-f b AQ(j,k) 



fa 



fb 



4tt 



(4) 



Here subscripts j and k distinguish each cell of the sky, 
Ail(j,k) denotes the solid angle of the (j,k) cell, and f a means 
the average of / calculated as 

fa = 2_^fa(j,k) — . (5) 



4tt 



By definition, 3 ranges from -1 to +1. When S = +1(— 1), the 
two distributions are exactly the same(opposite). When 3 = 0, 
we cannot find any resemblance between the two distributions. 
A source at a distance d, approximately contributes to arrival 
cosmic rays with the weight of d[~ 2 since all sources are as- 
sumed to have the same injection power. Therefore, f s (j,k) is 
calculated as ^ 1 jd? , where i runs over sources in the (j,k) 
cell. 

4. RESULTS 

For a given source distribution, we can simulate the arrival 
distribution of UHE protons using the calculation method ex- 
plained in the previous section. We investigate the correlation 
between the arrival distribution and the source distribution, and 
discuss the possibility to unveil the local structure by UHECRs. 

In FigQ] we demonstrate arrival distributions of UHE pro- 
tons above 10 19 8 eV predicted by a specific source distribution 
with 10~ 5 Mpc -3 in the galactic coordinate. The EGMF is not 
included. The source distribution within 100 Mpc is shown on 
its upper right panel. The radii of circles in this figure are in- 
versely proportional to source distances. The sources within 
50 Mpc are shown with bold circles. The number of simulated 
events are 50 (upper left), 200 (lower left), 500 (upper mid- 
dle), 1000 (lower middle), and 2000 (lower right). Note that 50 
events on all the sky correspond to the Auger result at present 
since Auger obse rves cosmic rays in the southern hemisphere 
dRoth etalJl2007h . 

When 200 events are detected, strong event clusterings from 
nearby sources (within 50 Mpc) can be observed. The distri- 
bution of more distant sources are not found yet except for the 
directions that sources are positionally concentrated. Detection 
of more than 500 events enables us to find event clusterings in 
the direction of almost all sources within 100 Mpc. Thus, we 
can find graphically that detection of about 500 cosmic rays 
above 10 19 8 eV makes us unveil nearby source distribution of 
them. 
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FIG. 1. — The arrival distributions of UHE protons above 10 19 8 eV predicted by a specific source distribution with 1(T 5 Mpc~ 3 (upper right) in the galactic 
coordinate. The EGMF is not included. The source distribution within 100 Mpc is shown as radii of circles inversely proportional to source distances. The sources 
within 50 Mpc are shown with bold circles. The simulated number of events are 50 (upper left), 200 (lower left), 500 (upper middle), 1000 (lower middle), and 2000 
(lower right). 



0J 
O 

o 
o 



S-h 
O 

o 



0.5 



0.5 



0.5 



TTj I I [ 1 1 T 1 1 1 I I I 1 1 1 I I [ llll| ITFTj I I I INI I I I I III T 

" (B,nj -(0.0,1 0- 4 ) " (R,n j-(0.1,10- 




~ (B,n s ) = (0.4,10- 4 ) 

x E>10 196 eV 
* E>10 198 eV 
■ E>10 W0 cV 




(B,n s )-(0.0,10- 5 ) " (B,n s )-(0.1,10- 5 ) " (B,n s )-(0.4, 1Q- 5 ) 




10 10 2 10 3 10 4 10 10 2 10 3 10 4 10 10 2 10 3 10 4 
Number of events 



FIG. 2. — Correlation coefficients between the arrival distribution of UHE protons above 10" 6 (crosses), 10 19,8 (triangles) and 10 20 °eV (squares), and their 
source distribution within 100 Mpc for several EGMF strengths and the number densities of UHECR sources as a function of expected number of cosmic rays on 
all the sky. The size of the cell for the calculation of the correlation coefficients is chosen to be 2° X 2°. The three marks on the same number of events are a 
little shifted horizontally for visibility. Units of B, normalization factor of the EGMF strength, and n s , the source number density, is fj,G and Mpc~ 3 respectively. 
The short (thick) error bars represent the fluctuations due to the finite number of observed events and the long (thin) error bars are total errors including the cosmic 
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We should discuss statistically (not graphically) the proba- 
bility that future observations can unveil UHECR source distri- 
bution. Figj2] shows predicted correlation coefficients between 
the arrival distribution of UHE protons and their source distri- 
bution within 100 Mpc for several EGMF strengths and source 
number densities as a function of expected number of cosmic 
ray events. The size of the cell for the calculation of the correla- 
tion coefficients is chosen to be 2° x 2° since angular resolution 
of UHECR experiments, 1°, is taken into account in simulating 
the arrival directions. The numbers of events that we calcu- 
late the correlation coefficients are 10, 20, 50, 100, 200, 500, 
1000, 2000, 5000, 10000 on all the sky for E > 10 19 6 (crosses), 
E > 10 198 (triangles), and E > 10 20 °eV (squares). The three 
marks on the same number of events are slightly shifted hori- 
zontally for visibility. The EGMF strengths, B, and the source 
number densities, n s , which are adopted for the simulation are 
represented on each panel. There are two errors. The thick 
(short) error bars show the statistical error, i.e., the fluctuations 
due to the finite number of observed events, averaged over 100 
realizations of the source distribution for a given source num- 
ber density. The random event selection to estimate the fluctu- 
ations is performed 100 times for each source distribution. On 
the other hand, the thin (long) error bars represent the total er- 
ror, i.e., the statistical error plus the cosmic variance. The latter 
error bars also include the variation between different realiza- 
tions of source distributions. Note that statistical errors on an 
observation are estimated as the thick error bars since only a 
source distribution is true in the Universe. 

First, we deliberate the correlation in the case of no magnetic 
field. A lot of features of the correlation coefficients common 
to finite EGMF cases are found in this simplest case. The three 
left panels show the correlation coefficients in the case of no 
EGMF. When the numbers of observed events increase, the cor- 
relation coefficients start to converge, and the final values can 
be estimated. The statistical errors (thick error bars) decrease. 
The cosmic variance remains at large number of events since 
it is uncertainty to originate from different source distributions. 
Note that a source distribution used in Fig Q] predicts correla- 
tion coefficients similar to the averages of those in the middle 
left panel of Figj2] The correlation coefficient cannot converge 
on 1.0 in spite of rectilinear propagation due to the finite an- 
gular resolution, the energy loss processes, and also because 
cosmic rays do not always come from sources within 100 Mpc. 
Since UHE protons from more distant source needs the higher 
energy at the source in order to reach the Earth with the same 
energies, such source contributes to arrival cosmic rays smaller 
than nearer source even if a factor of df 2 is taken into account. 
In other words, the arrival distribution of many cosmic rays is 
only approximately the same spatial pattern as their source dis- 
tribution weighted by a!, -2 . 

We focus on the behaviors of averages of the correlation co- 
efficients on each panel. The correlation coefficients with pro- 
tons above 10 20 eV predict larger values than those with lower 
energies at the small number of observed events. As the num- 
ber of detected events increases, the coefficients with highest 
energies converge on the smaller values than those with lower 
energies. The reason is GZK mechanism. UHE protons above 
10 20 eV can come from sources only within the GZK sphere 
(~ 50Mpc). On the other hand, protons with lower energies can 
arrive at the Earth outside 100 Mpc. Thus, the correlations are 
better at the small number of events. The difference between 
the correlation coefficients with protons above 4 x 10 19 eV and 



those above 10 eV is statistically significant at la level since 
strengths of errors on observation are estimated as the thick er- 
ror bars. Since the radius of the GZK sphere is less than 100 
Mpc, there are sources within 100 Mpc which cannot be traced 
by cosmic rays above 10 20 eV. Thus, as the number of observed 
events increases, the correlation coefficients with such highest 
energy protons converge on a little smaller values than those 
with lower energy protons. The correlations of protons with 
E > 10 19 ' 8 eV are better than those with E > 10 19 ' 6 eV, because 
lower energy protons can reach the Earth from more distant 
sources. 

Next, we investigate the correlation including EGMF. The 
three middle panels and the three right panels of Fig f2] show the 
correlation coefficients in the case of B = 0. 1 and 0.4^G respec- 
tively. The EGMF diffuses trajectories of UHEC Rs during their 
prop agation and obscures their arrival directions dTakami et al.l 
2006). Therefore, stronger EGMF predicts weaker correla- 
tion in any source number density. Compared to the cor- 
relation coefficient in the case of no EGMF, that of protons 
with the lower energy threshold is worse since lower energy 
protons are deflected by EGMF more largely. Protons with 
E > 10 196 , 10 19 8 eV predicts smaller correlation than those 
above 10 20 eV. 

The number of events that the correlation coefficients start 
to converge is weakly dependent on the EGMF strength and 
depends on the source number density. In the case of 
10~ 6 Mpc -3 , such numbers of events are about 1,000, 200, 
and 50 for E > 10 196 , 10 198 and 10 20 eV respectively. At 
larger number densities, such numbers increase. In the case 
of 10" 5 (10 -4 ) Mpc" 3 , these numbers are about 5,000(> 10000), 
500(5000), and 100(500) for£ > 10 19 6 , 10 198 and 10 20 eV re- 
spectively. Detection of such number of events enables us to 
unveil UHECR source distribution at visible level by UHE pro- 
tons, which concerned with the final value of the correlation 
coefficients. Both cosmic rays above 10 19 8 eV and those above 
10 20 eV are good indicators for unravelling the source distri- 
bution since their correlation values converge to nearly equal 
values for 10~ 4 and 10~ 5 Mpc -3 . However, the detected num- 
ber of cosmic rays above 10 19 8 eV reaches the number that the 
correlation coefficient starts to converge in shorter observing 
period than that of cosmic rays above 10 20 eV since cosmic rays 
above 10 20 eV are almost not detected in the presence of the 
GZK steepening. Note that our source model predicts the GZK 
steepening. Thus, UHE protons above 10 19 8 eV are best indi- 
cators to unveil the source distribution within 100 Mpc. Note 
that grou nd based detectors such as Auger and Telescope Ar- 
ray (TA) dFukushima et alj[2007f) can observe only about half 
hemisphere, so only half of such numbers are observed. 

Fig[3]is the simulation of the arrival distribution of UHE pro- 
tons above 10 19 8 eV in the case of B = 0.1 and 0.4/jG The 
source distribution is the same one in FigQ] In FigQ] and [3] 
the arrival distributions with 200 protons (less than 500 events, 
which is the number that the correlation coefficients start to 
converge) include several number of event clusterings. Com- 
pared to the three cases of the EGMF strength, the arrival direc- 
tions of protons are diffused by EGMF. We can also find event 
clusterings in the directions of nearby sources or source cluster- 
ings. Such event clusterings are not lost by a structured EGMF 
although their angular scales are spread. At 500 event detec- 
tion, the event clusterings increase and the local structure of the 
universe appears. 

The larger size of the cell for the calculation of the corre- 
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FIG. 3. — The arrival distributions of UHE protons above 10"' 8 eV predicted by the same source distribution used in Fig. 2. The upper panels show those taking 
into account propagation in the EGMF with B = O. l/iG and the lower panels are the same, but B = 0.4/iG. 
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FIG. 4.— Same as Fig|5] but for 4° X 4°. 
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lation coefficients is expected to lead to larger correlation co- 
efficients since EGMF obscures UHECR arrival directions as 
mentioned above. Fig|4] shows the same figure, but the size of 
the cell is chosen to be 4° x 4°. The correlation coefficients are 
larger than those in Figf2]on all panels. These increase even in 
the case of B = 0.0/iG since a part of arrival cosmic rays from 
a source is distributed in a small scale, but larger than 2° due 
to finite angular resolution of UHECR experiments. The values 
of the coefficients are slightly not increase even if the size of 
the cell is larger than 4° x 4°. In this figure, we find that the 
numbers of events starting the convergence of the correlation 
coefficients are almost unchanged. 

5. CONCLUSION 

In this paper, we calculated the arrival distribution of UHE 
protons taking into account energy losses and deflections by 
EGMF during propagation in intergalactic space in order to in- 
vestigate the possibility that future observations of UHECRs 
can unveil the local structure of UHE universe. In order to re- 
produce a realistic situation, we adopted a structured EGMF 
model and source distributions which reproduce the observed 
local structures. The arrival distribution of UHE protons was 
compared statistically to their source distribution using the cor- 
relation coefficients. As the number of observed events in- 
creases, the correlation coefficient increases and converges to 
some value which represents the ability to unveil the source 
distribution by UHE protons, i.e. charged particles. Thus, the 
number of events that the correlation coefficient starts to con- 
verge is an important number for UHECR observations. In 
other words, detection of such number of events allows us to 
unravel UHECR source distribution. We found that UHECRs 



above 10 19 8 eV are best indicators to decipher their source dis- 
tribution within 100 Mpc from discussion based on the final 
values of the correlation coefficients and GZK mechanism, and 
5000, 500, and 200 event detections above 10 19 8 eV on all the 
sky can unveil their source distribution for the source number 
densities of 10~ 4 , 10~ 5 , and 10~ 6 Mpc -3 respectively. Note that 
ground based detectors observe only about half hemisphere, so 
only half of such numbers are requested. 

In this study, we took only EGMF into account as magnetic 
field, i.e., neglected GMF GMF deflects trajectories of UHE 
protons efficiently by its regular components, which consist 
in spiral and dipole co mponents dAlvarez-Muniz etalJ. 12002; 
lYoshiguchi et alj l2003bh . A turbulent component of GMF 
very weakly change th e arrival directions of UHE protons 
dYoshiguchi et alj2004l) . The deflection angles of UHECR pro- 
tons are a few degree at around 10 20 eV except for the direction 
of Galactic Center. Such deflection disturbs the spatial pattern 
of UHECR arrival distribution at a few degree scale. The effect 
of GMF is one of our future investigations. 

A lot of inquiries on UHECR source number density result 
in around 10~ 5 Mpc -3 on ground of the AGASA results as in- 
troduced in section Q] If this number density is true, five year 
observation by Auger and future o bservation by TA and Ex- 
treme Universe Space Observatory dEbisuzaki et alj |2007) will 
reveal the distribution of nearby UHECR sources. The dawn of 
the UHE particle astronomy is just around the corner. 
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Scientific Research provided by the Ministry of Education, Sci- 
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